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Abstract 

We consider the massless two-loop two-point function with arbitrary powers of the propa- 
gators and derive a representation, from which we can obtain the Laurent expansion to any 
desired order in the dimensional regularization parameter £. As a side product, we show that 
in the Laurent expansion of the two-loop integral only rational numbers and multiple zeta 
values occur. Our method of calculation obtains the two-loop integral as a convolution prod- 
uct of two primitive one-loop integrals. We comment on the generalization of this product 
structure to higher loop integrals. 



1 Introduction 



Quantitative predictions from perturbation theory are crucially linked with our ability to calculate 
loop integrals. An object of extensive study has been the master two-loop two-point function, 
with massless internal propagators but arbitrary powers of the propagators. The name "the master 
two-loop two-point function" is justified, since all other two-loop two-point topologies can be 
obtained from this one. Allowing arbitrary powers of the propagators is important for three- 
or four- loop calculations, where the integration over self-energy insertions on the propagators 
shifts the power of the propagator from unity to 1 + e. (As usual we work within dimensional 
regularization with D = 4 — 2£.) 

In this paper we consider the massless two-loop two-point function. This integral is not 
only of practical importance from a phenomenological perspective, but received also quite some 
interest from number theorists. As far as the phenomenological side is concerned, this integral 
is implemented into the Mincer package dlU and used for example for the calculation of 
the total hadronic cross-section in electron-positron annihilation. From the number theoretical 
perspective there has always been the open question which types of (transcendental) numbers 
appear in the £-expansion of this integral. From explicit calculations it is known, that in the 
lowest orders multiple zeta values occur. 

The history for the calculation of this integral dates back, to the best of our knowledge, to 
1980, when the £°-coefficient for a specific (non-trivial) combination of powers of propagators 
was calculated using the Gegenbauer polynomial x-space technique |3|. In 1984 Kazakov ob- 
tained for all powers of the propagators of the form Vy = 1 + cij£ the result up to £ 3 , a year later 
the £ 4 -term followed 0|5l|6|. The momentum of pushing the £-expansion of this integral fur- 
ther was then taken up by Broadhurst and collaborators: In 1986 the £ 5 -term was calculated, in 
1988 followed the £ 6 term In the mean time the original Gegenbauer technique had been 

refined [9| and a representation of the integral for a specific combination of the powers of the 
propagators in terms of hypergeometric 3i*2-functions with unit arguments was found ll9l[T0l . The 
year 1996 brought an advance of two further terms in the expansion (£ 7 and £ 8 ) [ 10 1 . Finally, last 
year the £ 9 -term was announced ifTTl . The calculation of most of these terms relied on symmetry 
properties of the two-loop integral. It is known that this two-loop integral has the symmetry 
group Z2 x S(, llT2"l l71l8l. which is of order 1440, e.g. there are 1440 symmetry relations. The 
calculations cited above exploited the fact that the symmetry properties allow to reconstruct the 
result of the expansion up to order £ 9 for powers of the propagators of the form Vj = 1 + a ; £ from 
the result of the expansion, where two adjacent propagators occur with unit power. However, it 
is known that beyond this order the symmetry relations are not sufficient to determine the solu- 
tion |[10|. In view of the last point, the calculation of this two-loop integral cannot be considered 
to be solved to a satisfactory level. We quote from the latest publication of Broadhurst IfTTl : 

"It is one of the many scandals of our limited understanding of the analytical content of 
perturbative quantum field theory that, despite many years of intense effort, we still do not know 
whether multiple zeta values suffice for even the Taylor expansion of the two-loop integral." 

In this paper we proof that multiple zeta values are sufficient (theorem 14.11) . This is one of 
the main results of this paper and solves a long standing open problem. We also show how to 
obtain the £-expansion to arbitrary order for arbitrary powers of the propagators by calculating 
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this integral with the help of a new method. For an introduction to calculational techniques for 
multi-loop integrals we refer to H131ll4l[T51l . 

The second important result of this paper is the method we employed for the calculation of 
the two-loop integral: We obtain the two-loop integral as a convolution product of two primitive 
one-loop integrals. The convolution product can be evaluated in terms of nested sums 11611171 . 
and the £-expansion of the original integral is obtained from the £-expansion of these sums. 
Generalizations of this product structure can be useful for the extension of exisiting packages for 
two-point functions (like Mincer EIH) to four loops. In this paper we only briefly comment on 
this possibility and focus on the two-loop two-point function. In view of later applications for the 
calculations of higher orders we show that our method for the calculation can be implemented 
efficiently on a computer. 

This paper is organized as follows: In sect. E]we define the two-loop integral and summarize 
known facts about this integral. In sect. |3]we show that the two-loop integral can be written as 
a convolution product of two one-loop integrals. We also discuss the factorization for three-loop 
two-point functions. In sect. E]we use this product structure to obtain the e-expansion for the 
two-loop integral. In sect. |5]we report on the implementation of our formulae into a symbolic 
computer code. Sect, ^contains a summary and our conclusions. An appendix collects some 
useful formulae for integral transformations. 

2 Review of known results for the two-loop integral 

The object of investigation is the following five-propagator integral 

/ (2 ' 5) (m-e,Vi,V2,V 3 ,V4,V 5 ) = 

2\ v 12345 -2m+2e f d D k\ f d D k2 1 



-2 / 2\ v i2345-2m+2e f d K\ f 



i%D /2 j in D/2 (_ k 2y2 (_£|y 3 (_ fc2 y4 ^_^v 5 > 

(i) 

where kj, = £2 — P, £4 = k\ — p, k$ = &2 — k\ , D = 2m — 2e and 

r(l+e)r(l-e) 2 
Cr = r(l-2e) • (2) 

The superscripts for an integral indicate the number of loops I and the number of propaga- 
tors n. The prefactor in front of the integral is inserted for later convenience. It ensures that /( 2 ' 5 ) 
is independent of (— p 2 ) and avoids a proliferation of Euler's constant jg. The corresponding 
Feynman diagram for this two-loop integral is shown on the l.h.s of fig. [U It is for our purpose 
sufficient to assume that the exponents Vj are of the form 

Vj = nj + ajE, (3) 

where the nj are positive integers and the aj are non-negative real numbers. The £-dependence in 
the powers of the propagators arises from higher loop integrals by integrating out simple one-loop 
two-point insertions into the propagators of eq. 0. 
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Integration-by-part identities II181I191I20M relate integrals with different powers of the propa- 
gators. For example, from the triangle rule tlTHl we obtain the identities 



[(D-V235 -v 5 ) +v 2 2+ (1- -5-) +v 3 3+ (4- -5-)] 1^ 
[(D - v 145 - v 5 ) +vil+ (2- - 5-) +v 4 4+ (3- - 5~)] 7 (2 ' 5) 



0, 
0. 



(4) 



Here, the operators i + and i raise, respectively lower, the power of propagator / by one. The 
integration-by-part relations can be used to relate the integral where all propagators occur to a 
positive integer power to simpler topologies, e.g. where one propagator is eliminated. 
The integral in eq.© has the obvious symmetries 



(vi,v 2 ,v 3 ,v 4 ) 

(Vl,V 2 ,V 3 ,V 4 ) 



(v 2 ,vi,v 4 ,v 3 ), 
(v 4 ,v 3 ,v 2 ,vi). 



(5) 



However, there are more symmetries, which relate the integral to itself, up to prefactors of prod- 
ucts of Gamma functions. To discuss the symmetry properties of eq. ([TJ it is convenient to 
introduce the function F(vq, Vi, v 2 , V 3 , v 4 , V5) related to /( 2 ' 5 ) by 



/ (2 ' 5) (v ,vi,v 2 ,v 3 ,v 4 ,v 5 ) 
r(2v -3) 2 



where 



and 



(2vo-3)r(3-v ) 2 r(vo-l) 6 



v 6 



r(v y ) 



1/2 



7 7 (v ,vi,v 2 ,v 3 ,v 4 ,v 5 ; 



v 7 = 2v -v 235 , 

V 8 = 2V -Vi 45 , 



3Vq-Vi 234 5 



v 9 =v 345 -v , 

V (10) = v 125 ~ v 0- 



(6) 



(7) 



(8) 



Here and in the following we will use the short-hand notation like v, jy t = V,- + Vy +V* to denote 
sums of indices if {i,j,k} £ {1,2,3,4,5}. The function i 7 (vo,Vi,v 2 ,v 3 ,v 4 ,V5) is invariant under 
the symmetry group Zo x St. ll8l. The symmetric group is generated by the six-cycle 

(v ,vi,v 2 ,v 3 ,v 4 ,v 5 ) -»■ (v ,v 2 ,v 5 ,v 4 ,3vo-vi 2345 ,v 3 ) (9) 

and the transposition 

(vo,vi,v 2 ,v 3 ,v 4 ,v 5 ) -> (Vo,-Vo + Vi 4 5,V 2 ,-Vo+V 345 ,Vo-V 5 ,Vo-V 4 ). (10) 

The group Z 2 is generated by the reflection 

(v ,vi,v 2 ,v 3 ,v 4 ,v 5 ) -> (v ,vo-vi,vo-v 2 ,vo-v 3 ,vo-v 4 ,vo-v 5 ). (11) 

Note that in general the generators of the symmetry group do not conserve the positivity a 1 > 
in the parameterization of eq. ©. 
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3 The product structure of the two-loop integral 

In this section we show that the two-loop integral in eq. ^ can be written as a (convolution) 
product of two one-loop integrals. To this aim we define the following two one-loop integrals: 

d D k x 1 



t(\2)/ \ -ir 2\ v i4-" J + e f d k\ 1 



V4 1 



;«..3. (m _ e , V2 , V3 , V5W) = ^^^--70 ^^^ ,(12) 

where x = (— p 2 ) / (— £ 2 ), y = (—p 2 ) / (—#4) • Note that the integral /a 1 ' 3 ) depends on the kinematic 
variables — p 2 , —k\ and — k\ only through the dimensionless ratios x and y. The integral /t 1 ' 2 ) is 
easily computed as 

pafm-evi v 4 ) - c- 1 F( ~ m + £ + Vl4) F(m ~ £ ~ Vl )F(m ~ £ ~ V4) (13) 
[m £,v l5 v 4 ) - c r r(vi)r(v4) r(2m-2e-v 14 ) ' ° 3) 

The one-loop triangle integral 7( 1,3 ) can be written as a double Mellin-Barnes representation: 

Yi+/°° 72+!°° 

/ (1 ' 3) (m-£,V2,v 3 ,v 5 ;x,y) = -^—^ J do J Jxy"°^ T /( 1 ' 3 )(m-e,V2,v 3 ,v 5 ;x,o), 

Yj— i°° Y2 — 

(14) 

where the function /( 1,3 )(m — £,V2, V3, V5;x,o) is given by 
/( 1 ' 3 )(m-£,v 2 ,v 3 ,v 5 ;x,a) = c r ! 1 



r r(v 2 )r(v3)r(v5)r(2m-2e-v23 5 ) 

xT(-a)r(-a + m - £ - v 35 )r(-x)r(-x + m - £ - V25) 
xr(a + x-m + £+v 235 )r(a + x+v 5 ). (15) 

The integration contours are parallel to the imaginary axis, with indentations, if necessary, to 
separate the "UV"-poles (T(-o + ...), T(-x + ...), T(-a-x + ...)) from the "IR"-poles (T{o + 
...), r(x +...), r(a + x+ ...)). It should be noted that the function /^' 3 )(m-£,V2,v 3 ,V5;x,a) is 
the double Mellin transform in x and y of the original integral /( 1,3 )(m — £,V2,v 3 ,V5;x,y): 



00 00 



/ (1 ' 3) (m-£,v 2 ,v 3 ,v 5 ;x,a) = J dx j dy x x l y° 1 7 (1 ' 3) (m-£,v 2 ,v 3 ,v 5 ;x,y). (16) 



From eq. (fT3l and eq. (fl4l one obtains the two-loop integral as 

/( 2 ' 5 )(m-£,Vi,V2,V 3 ,V 4 ,V 5 ) = 

Yi+f°° j2+i°° 

j^—y; J do J Jx/ (1 ' 2) (m-£,vi -x,v 4 -o)/ (1 ' 3) (m-£,v 2 ,v 3 ,v 5 ;x,o). (17) 

Yi— i°° j2—i°° 




Figure 1: Factorization in terms of diagrams: The two-loop two-point function on the l.h.s. is 
equal to the insertion of a one-loop three-point function into a one-loop two-point function. The 
insertion occurs at the shaded vertex. 



In eq. ( fTTt the two-loop integral is obtained as a (double) convolution product of two one-loop 
integrals. This is pictorially shown in Fig. ([T]). At the level of Feynman diagrams the product 
structure is given as the insertion of a Feynman diagram into another Feynman diagram. It is 
a well known fact that convolution products can be turned into ordinary products by applying 
a suitable integral transformation. In the case at hand, eq. (fTTl) factorizes by performing two 
inverse Mellin transformations in Vi and V4. If one sets 

Y3+/00 Y4+'°° 

fi> 2 Xm-e,x,y) = J dVi J dV 4 x-^y- v 'P' 2 \m-e } v h v 4 ), (18) 

Y3— i°° Y4~ ,0 ° 

Y3+/00 Y4+'°° 

/( 2 > 5 )(m-e,V2,V3,v 5 ;x,y) = ^— ^ J dv x J dv 4 x^ l y^ A / (2 ' 5) (m-£,Vi,V2,v 3 ,v 4 ,v 5 ), 

Y3— i°° Y4 — ! °° 

one obtains the double inverse Mellin transform of the two-loop integral /( 2 5 ) as the product 
of the one-loop integral /( 1 ' 3 ) with the double inverse Mellin transform of the one-loop integral 

/(l, 2). 

/ (2 ' 5) (m-£,v 2 ,V3,v 5 ;x,y) = I {1 ' 2 \m-e,x,y) x / (1 ' 3) (m - £,v 2 , v 3 , v 5 ;x,y). (19) 

Eq. (flTt or eq. (fl9l) is the advertised factorization of the two-loop integral into two one-loop 
integrals. At the level of Feynman diagrams the product is similar to the insertion operation de- 
fined by Kreimer within the context of renormalization I2 2TI23II24I . In the context discussed here, 
insertions occur only at one specified place. This implies that the product is associative. Note 
that in general an insertion product, which allows insertions at several places, is not associative. 
The factorization property, e.g. that non-primitive graphs can be written as convolution products 
of primitive graphs, generalizes to higher loops. For the three-loop two-point functions there are 
three basic topologies, usually named the ladder ("LA") topology, the Benz ("BE") topology and 
the non-planar ("NO") topology. They are shown on the l.h.s of fig. ©. The ladder and Benz 
topologies are given as convolution products of three one-loop graphs: 
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Figure 2: Factorization at three loops: The first two topologies factorize into one-loop diagrams, 
whereas the last topology factorizes into a one-loop diagram and a two-loop diagram. A graph 
on the right side of the product operator * is inserted into the shaded vertex of the graph to the 
left of the product operator. 
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7 M' 8) ( m - e ' V l' V 2 5 V3,V4,V 5 ,V6,V7,V 8 ) 

/_ 2\Vl2345678-3 m +3e f d°k X f d°k 2 f d°k 3 * 1 

1 P ' J m D l 2 J m D l 2 J in / 2 }}, (-fc 2 )"' 



Yi+'°° Y2+'°° Y3+'°° Y4+«°° 

= -^^4 j doi J dxi J do 2 J Jx2/ (1 ' 2) (w-e,Vi-Xi,V6-Oi) 

X /( 1 ' 3 )(m-£,V2-X2,V 5 -02,V7;Xl,Oi)/ (1 ' 3) (m-£,V3,V4,V 8 ;X2,a 2 ). (20) 

The momenta are defined for the ladder topology by k$ = k 3 — p, k^ = k 2 — p, k§ = k\— p, 
h = k 2 — k\ and k% = k 3 — k 2 . For the Benz topology we have 

4£ 8) ( m - e > V l> V 2 5 V3,V4,V 5 ,V 6 ,V7,V 8 ) 

/_ 2\Vl2345678-3m+3£ /"A f d°k 2 f d°k 3 A 1 

1 P ] J in D l 2 J in D ' 2 J in D / 2 f = \ (^_ k 2^ 

Yi+i°° Y2+ 1 ' 00 Y3+'°° Y4+' 00 

j doi J d%\ J da 2 J Jx 2 / (1 ' 2) (m-£,vi -xi -x 2 ,v 5 -- 

Yi — i°° Y2 — i°° Y3 — ! °° Y4 — ! °° 

±1 1 

X 



= c r 3 (- 



-01) (21) 



X / (1 ' 3) (m-£,V2367 -m + £ + a2+X2,V4,V8-02;Xi,Oi)/ (1 ' 3) (m-£,V2,V7,V 6 ;X2,02). 

The momenta are defined for the Benz topology by k<\ = k 3 — p, k$ = k\ — p, k^ = k 2 — k\, 
k-j = k 3 — k 2 and k% = k 3 — k\ . The non-planar topology factorizes only into a one-loop two-point 
function and the crossed two-loop three-point function: 

4 3 ( j 8) (m-£,Vi,V 2 ,V3,V4,V 5 ,V 6 ,V7,V 8 ) 

-3 / 2^12345678-3m+3e f d°k X f d°k 2 f d°k 3 ^ 1 

" Cr [ P ] J m D l 2 J in°/ 2 J in D / 2 } = \ ( ] 

Yl+i'oo Y2+i°o 

= J do J ^ 7(1 ' 2) ( m -e 5 Vi-x,v 6 -o)4 2 ' 6) (m-£,V2,V3,V4,v 5 ,V7,v 8 ;x,o). 

Yi-ioo Y2-' 00 

Here, I x ' ; is the double Mellin transform 

4 2,6) (m-£,V2,V3,V4,V 5 ,V7,V 8 ;X,0) = 

OO CO 

j dxj dyx x ~ l y ~ l 4 2 ' 6) (m-£,V2,v 3 ,V4,v 5 ,V7,v 8 ;*,;y) (23) 



of the two-loop crossed vertex function 

4 2 ' 6) (m-£,V2,V 3 ,V4,V 5 ,V7,V 8 ;^v) = (24) 



-2 (_ 2\ v 234578-2m+2e f d k2 f d £3 1 

Cr 1 p ' J i% D ' 2 J i%°i 2 (-qy 2 {-k 2 3 y 3 (-k 2 4 y 4 (-k 2 5 y 5 (-k 2 ) Vi (-*§)* 

The momenta for eq. (l22b and eq. d24b are defined by k^ = k^ — p, k$ = k\ — ki + k^ — p, 
k(, = k\ — p, kj = k2 — k\ and ^8 = ^3 — &2- The variables x and y are defined by x = (—p 2 ) / (—kf) 
and y = (—p 2 ) / (— &|). For all topologies the factorization is shown pictorially in fig. ©. 



4 Evaluation of the two-loop integral 

Eq. (fT7t is the starting point for the further evaluation of the two-loop integral: 



Yi+i°° Y2+' 00 

1 f f r(-a)r(-o + m-E- v 35 )r(o + m-£-V4) 



/(2,5) =c ^ dO dX 



(2m) 2 J J r(-o+v 4 ) 

Yi — i'c<3 Y2 — '°° 

r(-T)r(-t+m-e— V25)r(t + m-e — Vi) 

x r(-x+vi) 

r(— a - x — m + £ + Vi 4 )r(o + x - m + £ + V235)r(a + x + v 5 ) 
r(o + X + 2m-2£-Vi 4 ) 



(25) 



with 



c = r(v 2 )r(v 3 )r(v 5 )r(2m - 2e - v 235 ) ' (26) 

The strategy for the evaluation is as follows: We first close the contours and evaluate the integrals 
with the help of the residuum theorem. This technique has a long history, a recent example is the 
calculation of the triple box ll25l . The semi-circles at infinity needed to close the contours give a 
vanishing contribution provided that 

vi +V125 -2m+2e < 1, 
V4+V345 -2m + 2£ < 1, 

-1 < (vi+Vi25-2m + 2e) + (v 4 +V345-2m + 2e). (27) 

From the residues we immediately obtain nested sums. All these sums can be brought to a stan- 
dard form and then expanded into a Laurent series in £ with the help of the algorithms described 
in ifTol . We use the program "nestedsums" ifTTI and the "GiNaC'-library [26 1 for this purpose. 
The algorithms A and B of lfT6l are generalizations of algorithms for the manipulations of har- 
monic sums and harmonic poly logarithms, described in H27U28ll2"9l . 

Closing the contour of the a-integration to the right, one picks up the residues of T(— o), 
r(-o + m- £ — V35) and r(— o — x — m + £ + Vi4). All other a-dependent Gamma functions in 
the numerator have poles to the left of the contour and therefore do not contribute. The basic 
formula for the residuum of Euler's Gamma-function reads: 

(-1)" 

res (T(-x + a),x = a + n) = -- — J— (28) 
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The two-loop integral can therefore be written as a sum of three terms, 

/(2,5) = Til) + Ti2) + T(?)i (29) 

where each term is obtained by taking the residues of one Gamma function from the set T(— o), 
r(-G + m-£-V35) and r(-G-X-m + £ + Vi4). Explicitly, 

w n to n\ 2%i J r(-x+Vi) W ' 



n=0 

J2-1 



with 



r(-n + m-£-V35)r(n + m-£- V4) 
r(-n+v 4 ) 

r(— x — n — m + £ + Vu)r(x +n — m + £ + V235)T(x + n + V5) 



x 



H (2) 



r(x + n + 2m-2£-Vi 4 ) 
r(— n - m + £ + V 35 )r(n + 2m - 2£ - V345) 

r(-n-m + £ + V 3 45) 
r(— T - n - 2m + 2£ + Vi3 45 )r(x + n + V2)T(x + n + m - £ - V3) 



x 



*<3) 



T(x + n + 3m - 3£ - V1345) 
r(n - 2m + 2£ + Vi 2 345)r(rc -m + £ + V145) 
T(n + m — e) 

r(x-n + m-£-Vi4)r(x-n + 2m-2£-Vi34 5 )r(-x + n + Vi) 

x — — . w-U 

r(x-n + m-£- Vi) 

The procedure is then repeated by closing the x-integration contour to the right. For example, 
for the term 7m we have to evaluate the residues of T(— x), T(— x + m — £ — V25 ) and T(— x — 
n — m + E + Vu). The residues of r(x + m — £ — Vi), r(x + n — m + £ + V235) or r(x + n-fVs) are 
always to the left of the contour. In summary, the two-loop integral can be written as the sum of 
the following terms: 

/( 2 > 5 ) = r (1;1 ) +7(1,2) +7(1,3) + 7(2,1) + T {2,2) + 7(2,3) + 7(3,1) + 7(3,2) + 7(3,3) + 7(3,4) + 7(3,5)- 

(32) 

Table © shows the correspondence between each term Tr^n and the Gamma functions, from 
which the residues are taken. The term 7(3,3) corresponding to the residues in T(— x + n + Vi) 
gives no contribution. This Gamma function is always accompanied by 1 /T(— x + Vi) . Since 

n r 7i!tT' = (-, + v I+ „-l)...(-x + v 1 ) (33) 

r(-x+vi) 

is free of poles we have 

7(3,3) = 0. (34) 



10 





a 






X 






r( 


-a) 




r( 


/ 


r (l,2) 


r( 


-a) 




r( 


-X + m-£- V25) 


r (l,3) 


r( 


-a) 




n 


— X — n — m + £ + Vi 4 ) 


To i i 


r( 


— G + m - 


-£ — V35) 


r( 


/ 


7}? 91 


r( 


— a + m - 


-£- V35) 


r( 


-X + m-£-V25) 




r( 


— a + m - 


-£ — V35) 


r( 


— X — n — 2m + 2e + V1345) 


r (3,l) 


r( 


— O — T - 


m + £ + Vi 4 ) 


r( 


-x) 


7(3,2) 


r( 


— O — T - 


m + £ + Vi 4 ) 


r( 


-X + m-£-V25) 


r (3,3) 


r( 


— O — T - 


m + £ + Vi 4 ) 


r( 


-X + n+Vi) 


r (3,4) 


r( 


— O — T - 


m + £ + Vi 4 ) 


r( 


X — n + m — £ — Vi 4 ) 


r (3,5) 


r( 


— O — T - 


m + £ + Vi 4 ) 


r( 


X - n + 2m - 2e - V1345) 



Table 1: Correspondence between the terms T^ k j^ and the Gamma functions, from which the 
residues are taken. 



Furthermore, we can show that the terms 3) and To, 4) as well as the terms 7^3) and To, ,5) 
cancel each other: 

7(i,3) + 7(3,4) = 0, 7(2 5 3) + r ( 3 55) = 0. (35) 

To present the results after all residues have been taken, we introduce two functions G± with ten 
arguments each: 

G±(a]_,a 2 ,a 3 ,a4;b]_,b 2 ,b 3 ;C]_,C 2 ,C 3 ) = 

f y (-ir +j r(Tn-j-a l )r(±n + j + a 2 )r(±n + j + a 3 ) 

^ T(TnTbi)T(n + b2)n-j-c l )T(j + c 2 ) 

T( T n T b 3 ) r(-j-c 3 ) ' (M) 

together with two operators L d and acting on the arguments as follows: 

L d G± (a h a 2 , a 3 ,a 4 ; b\ , b 2 , b 3 ; c\ , c 2 , c 3 ) = 

G±{a\ + d,a 2 + d,a 3 + d,a 4 + d;bi + 2d,b 2 + d,b 3 + d;c\,c 2 ,c 3 ), 
%dG± {ai,a 2 ,a 3 ,a 4 \b\MM\c\,c 2 ,c 3 ) = 

G±(oi + d, a 2 + d,a 3 + d, a 4 + d; b\ , b 2 , b 3 ; c\ + 2d, c 2 + d,c 3 +d). (37) 

Then 

/(2,5) _ c (l_|_^ m _ e _ V35 +^„_ e _ V25 + X ra - e -v35^m-e-V25) 
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G + (m - £ — V14, —m + £ + V235 , V5, 2m - 26 — V14; —m + £ + V35 , m - £ - V4, -V4; 

-m + £ + V25,m-£- Vi, — Vi) 

+c(l + ^ M _ e _v 25 )G_(-Vi,m-e-Vi4,2m-28-Vi3 45 ,m-8-Vi; 

-2m + 2£ + V12345, — m + £ + V145, m - £; — m + £ + V25, m - £ - Vi, -Vi) . (38) 

Here, 1 denotes the identity operator with a trivial action on the arguments of the functions G±. 
To proceed further, we now show how to transform the functions G± to a standard form, such 
that they can be expanded in £ with algorithms of ifTSlfTVl . Since the shift operators L and % 
modify only the arguments but not the structure of the functions G±, it is sufficient to discuss one 
example for each function. For the function G + we discuss as an example the first term (without 
any application of the shift operators L or $Q of eq. d3~5b . which corresponds to the term Tn^y 

r (1;1) = cG + (m-£-Vi4,-m + £ + V235,V5,2m-2£-Vi4;-m + £ + V35,m-£-V4,-V 4 ; 
-m + £+V25 ; m-£-V],-Vi) 
" y (_i)("+j)r(-n-j-m + £ + Vi 4 )r(n + j-m + £ + V235)r(n + j+v 5 ) 
C ~oh nljl r(n + j + 2m-2e-Vi 4 ) 

r(-n + m-£-V3 5 )r(n + m-£-v 4 ) T(-j + m-e- V 2 5)r(j + m-£-Vi) 
X T(-n+v 4 ) n-j+Vi) ^ } 

If Vi and V4 are not positive integers, one may use the reflection formula 

r(-n+x) = r(x)r(i-x) ( ~ 1)w (40) 

1 yn + 1 — x) 

to flip the Gamma functions, where a summation index occurs with a negative sign. One obtains 
in that case 

r (M) = cr(-m + e+Vi4)r(l+m-e-Vi 4 ) 

r(m-£- V 35 )r(l -m + £ + V 35 )r(m-£- V25)r(l -m + £+V25) 

x r(vi)r(i-vi)r(v 4 )r(i-v 4 ) 

x y y r(n + m-£-V4)r(n+l-v 4 ) T(j + m-E-Vi)T(j + 1 - Vi) 
X hh r(n + l)r(/» + 1 - m + e + v 35 ) T(y + l)r(y + 1 - hi + e + V25) 
r(n + j'-m + £ + v 2 35)r(n + 7+V5) 



(41) 



T(n + j' + 2m-2£-Vi 4 )r(n + j + l+m-e-Vu) 

This is a double infinite sum with unit arguments and can be evaluated with algorithm B of lfl6l 
im. The sum in eq. d4*TT) is a generalization of the first Appell function. The Laurent expansion 
to any fixed order will contain only rational numbers and multiple zeta values. If Vi or V4 are 
positive integers, the infinite sums over j or n terminate. For example, if Vi is a positive integer, 
while V4 is not, one obtains 

r(m-£-V 35 )r(l - OT + £ + V 35 )r(-m + £ + Vi4)r(l + m-£-Vi4) 

(M) = c r(v 4 )r(i-v 4 ) 



12 



v ^ r(-j + m-e-v 2 5)rX/ + m-£-Vi) " r(n + m-£-v 4 )r(n + l -y 4 ) 

x Po r(j + i)r(-j+vO ^ r(n + i)r(n + i-m+£+v 35 ) 

x T(n + j-m + £ + V 2 3 5 )r(n + j + V 5 ) 

X T(n + ;' + 2m-2e-Vi4)r(n + 7 + 1 +m-e- Vw) ' 

This is a finite sum of single infinite sums, more precisely, it is a sum of Vi terms, each containing 
a 4F3 hypergeometric function with unit argument. These are evaluated with algorithm A of 
ifToirm . Again, the Laurent expansion to any fixed order will contain only rational numbers 
and multiple zeta values. The case where V4 is a positive integer while Vi is not, is completely 
analog. If both Vi and V4 are positive integers, one obtains 

VV (-1) ( " +J) r(-n-j+Vi4-m + £)r(n + j-m + £ + v 2 35)r(n + j+v 5 ) 
(M) °hh nljl r(n + j + 2m-2e-v u ) 

r(-n + m-£-V35)r(n + m-£-V 4 ) r(-j + m-£-V2 5 )r(j + m-£-Vi) 
X T(-n + v 4 ) r(-y+Vi) ' 1 } 

e.g. a finite sum of products of Gamma functions. This is rather trivial and the Laurent expansion 
contains only rational numbers and zeta values. 

We now turn to the function G . As an example we discuss 

r (3il ) = cG (-Vi,m-e-Vi4,2w-2e-Vi345,m-e-Vi; 

-2m + 2e+v 12345, — m + e+Vi45,m — e;— m + e+V"25,m — e — vl,— Vi) 
r(n-2m+2e+Vi2345)r(n-m + 6+Vi4 5 ) 



00 00 



cY V 

tof^o nljl r(n + m-e) 

r(-j + m-£-v 25 )r(j + m-£-Vi) 

X — 



r(-j+vi) 

r(-n + j + m-£-Vi 4 )r(-n + j + 2m-2£-Vi345)r(n-7+Vi) 
r(— n + j+m — e— Vi) 



(44) 



If Vi is not a positive integer, we split the summation region into the regions j < n and n < j and 
subtract the double counted diagonal j = n as in 

00 00 co fi 00 J 00 

EE*W) = EE F (*>;)+EE F (*>-/')-E FM - (45) 

n=0 j=0 n=0 j=0 ;'=0n=0 «=0 

We obtain 

7( 3a ) = cT(m- 8 -Vi 4 )r(l -m + £ + Vi 4 ) 

r(m-e-V25)r(l -m + £ + v 2 5)r(2m-2£-Vi3 45 )r(l -2m + 28 + Vi3 45 ) 

x r(vi)r(i -vi)r(m-£-vi)r(i -m+£+vi) 

x y y r(n + Vi)r(n + l-m + £ + Vi) T(j + m-£-v l )T(j + 1 - vQ 

X „to E Q r(n + 1 - m + £ + V M )r(n + 1 - 2m + 2£ + V1345) T{j + l)T(j + 1 - m + £ + V 25 ) 
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r(w + j-m + e + Vi45)r(yz + j-2m + 2e + Vi2345) 
F(n + j+l)T(n + j + m-E) 

, v Wl , , v f f- r(/i + m-£-Vi 4 )r(n + 2m-2£-Vi34 5 ) 

+cT m-£-v 25 r l-m + £ + v 25 L L r / xi Z r i ^777 

„tojto r(n+l-Vi)r(n+m-e-Vi) 

r(j-m + 8+Vi4 5 )r(7-2m + 28 + V 12 345)r(n + 7 + m-8-Vi)r(n + 7 + l-Vi) 



x 



— c 



r(j + l)r(j+m-e) r(n + 7 + l)r(n + j + l-m + e+v 25 ) 

T(m - £ - v 25 )r(l -m + £ + v 25 )r(m - £ - Vi 4 )r(2m - 2£ - V1345) 

r(i-vi)r(m-e-vi) 

y, r(n + m-£-Vi)r(n-m + £+v 1 45)r(n-2m + 2£ + Vi 2 34 5 )r(n+l -Vi) 
X to r(n+l)r(n + l)r(n + m-£)r(n+l-m + £ + v 25 ) <0> 

The sums are of the same type as discussed in eq. (|4TTl and eq. (|42l . therefore the Laurent 
expansion contains only rational numbers and multiple zeta values. 

If Vi is a positive integer, we observe that in eq. (03) we have the building block 

r r(-j + vi) = ( v l-i')(Vl-J + l)-(Vi-J + n-l). (47) 

Therefore we obtain a non-vanishing contribution only if this sequence does not contain zero. 
This is the case for 

Region I: ;*<Vi, 

Region H: j > n + V\ - 1. (48) 

Region I contains only a finite number of residues from T(— x), whereas region II contains an 
infinite sum. We obtain 

r (3)1) =c(-l) v T(m-£-v 125 )r(l-m + £ + Vi 25 ) 

r(n + m — £ - V4)T(n + 2m-2e- v 345 ) T(j -m + e + Vu 5 )T(j - 2m + 2z + Vi 23 4 5 ) 



x 



x 



00 00 



II 



to ; to r(n + l)T(n+m- e) T(j + l)T(j + m - e) 

F(n + j+l)T(n + j + m-E) 
F(n + j'+l+Vi)r(n + ; + l-m + e + Vi 25 ) 
r(m-e-Vi 4 )r(l -m + e+Vi4)r(2m-2e-Vi345)r(l-2m+2e+Vi34 5 ) 
+ ° r(m-£-Vi)r(l-m + £+Vi) 

r(-j + m-£-v 25 )r(j + m-£-Vi) y T{n-m + £ + Vi 45 )r(n - 2m + 2e + Vi 234 5) 

x r(y+i)r(-7+vi) „t* r(n + l)r(n+m-e) 

x r(/z-j+vi)r(/i-j + i-m+e+vi) 

r(n- 7 + 1 -m + £ + Vi 4 )r(n- 7 + 1 -2m + 2e + V1345) ' 

Once again, the sums are of the same type as discussed in eq. (|4*TT) and eq. (|4~2"1) . and the Laurent 
expansion contains only rational numbers and multiple zeta values. 
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To summarize, we are able to write the two-loop integral /( 2 ' 5 ) as a combination of terms, 
which can be expanded in e to arbitrary order with the help of algorithms A and B of H161I17I . 
All these terms occur with unit arguments, therefore the Laurent expansion of /( 2 ' 5 ) involves only 
rational numbers and multiple zeta values. In particular, in the representation derived here, no 
alternating Euler sums occur. This answers a question raised recently by Broadhurst ifTTl : 

Theorem 4.1 Multiple zeta values are sufficient for the Laurent expansion of the two-loop inte- 
gral ft 2, 5 \m — £, Vi, V2, V3, V4, V5), if all powers of the propagators are of the form Vj = nj + afi, 
where the nj are positive integers and the aj are non-negative real numbers. 

We have shown above, that the Laurent expansion of the two-loop integral /( 2 > 5 ) can be expressed 
in multiple zeta values if the inequalities eq. dT7t are satisfied. These inequalities ensure that the 
semi-circles at infinity give a vanishing contribution. If for a specific combination of powers 
of propagators these inequalities are not satisfied, one may use the integration-by-part identities 
eq. © and the symmetry relations eq. © and eq. © - dTTT) to express the original integral as 
a linear combination of integrals, which fullfill the inequalities. The coefficients of this linear 
combination are either rational numbers or Gamma functions, which in turn expand into zeta 
values. Therefore the theorem holds also in the general case. 

The restriction to aj > in the parameterization v 7 = n; +ajE ensures that the different 
cases we treated above (e.g. whether the argument of a Gamma function is an integer or not) 
are sufficient. This condition could be relaxed at the expense of a more extensive case study. 
However, in all practical calculations the integrals occur with aj > 0. 

Finally we remark on a technical detail: Our method of calculation expresses the two-loop 
integral as a combination of several sums with unit arguments. Although the final result is finite, 
it is not guaranteed that all individual sums are convergent. Furthermore, the algorithms used 
for the Laurent expansion rely on partial fractioning. This may split a convergent sum into two 
divergent pieces, as illustrated by the example of the convergent sum 

00 1 

I^-Zn = L (50) 

„=i n{n+l) 

Partial fractioning splits 

1 1 1 



ran u n (51) 

n{n + l) n n+l 

The problem is easily circumvented by introducing a finite upper summation limit N or by mul- 
tiplying the summand by x n and taking the limit x — > 1 in the end. With both methods, the 
divergent pieces cancel at the end of the day and one obtains the correct and finite result. Within 
the second method, one inserts 

x a x x x 2 (52) 

into the integrand of eq. d25l) . Then all sums are convergent, provided < x\ < 1, < xi < 
1 and X2 < x\. The Laurent expansion for this more general expression will contain multiple 
polylogarithms in x\ and x^. The result for the two-loop integral /( 2 > 5 ) is recovered by first taking 
the limit x\—*\ and then the limit X2 — > 1. Note that the order of the limits cannot be exchanged. 
It is easily seen that in the limits the multiple polylogarithms reduce to multiple zeta values. 
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0(E) 





1 


2 


3 


4 


5 


6 


7 


weight 


3 


4 


5 


6 


7 


8 


9 


10 


time 


29 


54 


99 


185 


375 


910 


2997 


11741 


memory 


6 


7 


8 


12 


30 


104 


397 


1970 



Table 2: CP£/ rime in seconds and required memory in MB for the expansion o//( 2 ' 5 ) (2 — £, 1 + 
e, 1 + £, 1 +£, 1 + £, 1 +e) up to the indicated order/weight on a PC (1.6 GHz Athlon with 2 GB 
RAM). 

5 Results, checks and performance 

We have implemented the results of the previous section into a C++ program, which calculates 
symbolically the Laurent expansion of the two-loop integral for a user-specified set of parameters 
(m,Vi, ...,V5) up to the desired order in £. The program uses the "nestedsums"-library ifTTll and 
the "GiNaC-library ll26l . To simplify our results we use the Grobner basis for multiple zeta 
values provided in ll30l . 

To check our implementation we have written two independent programs and verified that 
they agree. Further we have checked that for trivial cases (e.g. all powers of the propagators 
integers) we obtain the correct (known) result. In addition, we have verified the symmetry rela- 
tions eq.©. Finally we compared for non-integer powers of the propagators with known results 
from the literature. As an example we quote the result for one particular integral: 

(1 - 2e) / (2 ' 5) (2 - £, 1 + £, 1 + £, 1 + £, 1 + £, 1 + £) = 6^3 + % A z + 372^ 5 £ 2 

+ (915^ - 864C|) £ 3 + (18450^7 - 2592^3) e 4 + (50259^ - 76680^ 3 - 2592^ 2 ) £ 5 
+ (905368^9 - 200340^ 6 ^3 - 130572^ 4 + 66384^) £ 6 

+ (2955330^io - 68688^,2 - 3659904^^3 - 1777680(g + 298728^) e 7 + O (£ 8 ) . 

(53) 

For multiple zeta values we use the notation 

oo j 1 

m\,mi = 2^ ~mY~m2' (54) 
«!>n 2 >0 n l H 2 

We have compared our result up to weight 10 with ifTTJl and found agreement. The order to 
which we can expand the two-loop integral in the parameter £ is not limited by our method of 
calculation. The only restriction arises from hardware constraints (e.g. available memory and 
CPU time). Table © shows for the case of the integral /( 2 ' 5 ) (2 - £, 1 + £, 1 + £, 1 +£, 1 + £, 1 + e) 
the dependence of the required CPU time and the required memory on the order of £ to which the 
Laurent expansion is calculated. In table © we also indicate the highest weight of multiple zeta 
values, which occur within a given order. Note that the weight is a more accurate measure of the 
complexity of a calculation, since in the example discussed here individual terms start at 1/e , 
but in the sum the coefficients of the terms l/£ 3 , l/£ 2 and 1/e vanish. In general we expect terms 
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up to weight 21 — 1 to occur in the finite part of a /-loop two-point function. We see from table 
that up to weight 8 our implementation is rather efficient. For higher weights the main limitation 
is given by the available memory. However, as far as practical applications are concerned, weight 
7 is sufficient to extend the exisiting (a 3 ) calculation H31 11321 for e + e~ — > hadrons to order a 4 . 
Here one would need the following integrals: 

(1-28) 7 (2 ' 5) (2 - e, 1+ 2e, 1 , 1 , 1 , 1) = 6& + % 4 £ + 142^ 5 £ 2 + (340^ 6 - 158^1) e 3 
+ (3034^7 - 474^3) e 4 + f^Cs - 6172^ e 5 

+ 



V 4 

& _ 9258^ 4 - 14640^3 + e 6 + (e 7 ), 



( 1 - 2e) / (2 ' 5) (2 - e, 1, 1, 1 , 1 , 1+ 2e) = 6^3 + % 4 e + 192^ 5 e 2 + (465^ 6 - 168g) £ 3 

+ (4509^ - 504^ 3 ) e 4 + (^f^s - 1620^,2 - 3252^) e 5 

+ (98490^9 - 14598^^4 - 15390^3 + 2676^) £ 6 + o (e 7 ) , 

(1-28) 7 (2 ' 5) (2 - 8, 1 + 8, 1 + 8, 1, 1, 1) = 6^3 + 9^ 4 e + 132^ 5 e 2 + (3 15^ 6 - 144g) e 3 

+ (2634^7 - 432^3 ) e 4 + (7749^8 - 5256^ 3 ) e 5 

+ (53160^9 - 12420^ 6 ^3 - 7884^ 4 + 1872g) e 6 + o (£ 7 ), 

(1-28) P' 5 ) (2 - 8, 1 + 8, 1 , 1 + 8, 1 , 1) = 6^3 + 9^ 4 e + 1 27^ 5 £ 2 + - 173^ 

+ - 519C 4 C 3 ) e 4 + - f - 5839 C5 C 3 ) s 5 

+ - 14340«s - + ^S) + O (£ 7 ), 

(1-28) 7 (2 ' 5) (2 - 8, 1 + 8, 1 , 1 , 1 + 8, 1) = 6^3 + 9^ 4 e + 132^ 5 e 2 + (3 15^ 6 - 204g) £ 3 
+ (2634^7 - 612^3) e 4 + f^^g - 7476^1 e 5 



£ 3 



2 

+ (53160^-17670^3 - 11214^^ + 3612^) £ 6 + o(£ 7 ), 

(1-28) (2 - 8, 1 + 8, 1 , 1 , 1 , 1 + 8) = 6£ 3 + 9^ 4 e + 1 57^ 5 £ 2 + - 179g) £ 3 

/26657>. \ 4 /124899>, 1215„ , \ 5 

+ ( — g— b - 537^^3 ) e 4 + [—^- ^ ~ ~2~^ 2 ~ 5521 CsC3 J e 5 

/1657525,. 1CA , cr r 23853.. 8776 r3 \ 6 , 7 , 
+ ( 24 ^9 - 15945C 6 C3 - -^—^5^4 + — Q J e 6 + (£ 7 ). (55) 

These integrals are in principle only needed up to order £ 4 for the a 4 corrections. However, 
knowing these integrals to higher orders can simplify integral reduction algorithms. For example, 
one may allow spurious poles in integration-by-parts identities. The CPU time and the required 
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memory needed to calculate each of these integrals is slightly below the corresponding numbers 
indicated in table ©. 



6 Summary and conclusions 

In this paper we calculated the £-expansion of the massless master two-loop two-point function 
with arbitrary powers of the propagators. We showed that to all orders in £ rational numbers 
and multiple zeta values are sufficient to express the result. Our method of calculation obtained 
the two-loop integral from a convolution of two one-loop integrals. We also discussed the cor- 
responding factorization for three-loop two-point functions. Finally we demonstrated that our 
method can be implemented efficiently on a computer. 
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A Integral transformations 

In this appendix we briefly summarize the Laplace, Mellin and Fourier integral transformations. 
Let fit) be a function which is bounded by an exponential function for t — ► ±00, e.g. 



1/(0 1 < Ke Cot for t 
\f(t)\<K'e- Clt fort 



00 



— 00, 



(56) 



Then the (double-sided) Laplace transform is defined for cq < Re o < c\ by 




(57) 



The inverse Laplace transform is given by 



Y+i'oo 




(58) 



The integration contour is parallel to the imaginary axis and Co < Re y < c\ . 
Let h(x) be a function which is bounded by a power law for x — > and x — > ■ 



00, e.g. 



\h(x)\<Kx- c ° forx^O, 
\h(x)\ <K'x Cl for x^ 00 



(59) 
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Then the Mellin transform is defined for cq < Re a < c\ by 

oo 

h M (a) = fdxh(x)x a - 1 . (60) 



o 

The inverse Mellin transform is given by 



h(x) = — J doh M (o)x °. 



(61) 



y-i° 



The integration contour is parallel to the imaginary axis and cq < Re y < c\ . There is a close rela- 
tion between the Laplace and the Mellin transform: If fit) = /z(e~ f ), then f L (a) = h u (c). Note 
that the inversion formulas can be obtained from the known inversion formula for the Fourier 
transform. To a function fit) we associate the Fourier transform 

oo 

f T (u) = J dtf(t)e- 2niut . (62) 

— oo 

The inverse transform is given by 

oo 

f{t) = [ duf r iu)e 2niut . (63) 



The Laplace transform and the Fourier transform of a function fit) are related by 

f ^ = *(si)- <64) 

As an example for the Mellin transform we consider the function 

h(x) = . X . (65) 
with Mellin transform h M (a) = T(— a)T(a + c)/r(c). ForRe(— c) <Rey< Owe have 

* c i /■ ^r(-o)r(o+c) _ a 



— / do — — ' \ y (66) 



(i+jc) c 2to 7 r(c) 

y— /'oo 

From eq. (l66t one obtains with jc = 5/A the Mellin-Barnes formula 

y+i'oo 



{A+B ^ = 25 7 rfg A B~ . (67) 



y— zoo 
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We often deal with integrals of the form 



y+foo 

im J r(o+c 2 )...r(o+c p ) r(-o+di)...r(-o+d 9 ) 

y— ;'oo 

If max (Re ( —a i ),..., Re( —a m )) < min (Re (&i ),..., Re(fc„) ) the contour can be chosen as a straight 
line parallel to the imaginary axis with 

max(Re(-ai),...,Re(— a m )) < Rey < min(Re(&i), ...,Re(fr„)) , (69) 

otherwise the contour is indented, such that the residues of T(a + ai), T{o + a m ) are to the 
right of the contour, whereas the residues of T(—a + bi), T(—a + b n ) are to the left of the 
contour. We further set 

a = m + n — p — q, 
P = m — n — p + q, 

(in n p q \ 1 

L «.i ' L h .i L O L (l ! )-~(m + n-p-q). (70) 
j=l j=\ j=l j=\ J L 

Then the integral eq. d68t converges absolutely for a > lf33l and defines an analytic function in 

|argx| < min^7i,a^j. (71) 

The integral eq. d68t is most conveniently evaluated with the help of the residuum theorem by 
closing the contour to the left or to the right. Therefore we need to know under which condi- 
tions the semi-circle at inifinty used to close the contour gives a vanishing contribution. This is 
obviously the case for |jc| < 1 if we close the contour to the left, and for |x| > 1, if we close the 
contour to the right. The case |jc| = 1 deserves some special attention. One can show that in the 
case P = the semi-circle gives a vanishing contribution, provided 

I < -1. (72) 

To derive this result, Barnes asymptotic expansion of the Gamma function for large x is useful: 

lnT(x + c) ~ ( x + c )lnx-x-W^-t ? 7^(--)\ (73) 

2 2n ^ n(n+l) \ x) 

where the Bernoulli polynomials B n (x) are given in terms of the Bernoulli numbers Bj by 

B n (x) = L("W"~ j - (74) 

7=0 V J / 
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